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1. Introduction 



At an early stage the universe was very likely a quark-gluon plasma. In heavy ion colliders 
we seek to recreate this state of matter and study its properties. Lattice gauge theory is ideally 
suited for the fully nonperturbative study of quantum chromodynamics under conditions close to 
thermal equilibrium. The insights gained from lattice simulations can be extrapolated through 
hydrodynamic modeling to the quasi-equilibrium expansion of the plasma. In this way lattice 
calculations provide crucial assistance in the interpretation of experimental results [||]. 

The baryon density was essentially zero in the early universe. At high densities, even at low 
temperature, popular tradition predicts that hadronic matter is also in a deconfined plasma state. 
At still higher densities even more unusual phases, such as a color-flavor locked phase have been 
proposed [Q]. For technical reasons, such high densities are beyond the reach of lattice simulations 
using standard methods. At least we can hope to simulate matter at the low densities found in heavy 
ion collisions. 

In this talk I consider only progress in calculations at zero baryon number density. In the 
companion talk Shinji Ejiri describes developments in nonzero density calculations [|J]. 

Here are highlights of recent advances at zero baryon density, which I will cover in this talk: 



• HotQCD study. A high statistics study is being carried out on the IBM BlueGene at the 
Lawrence Livermore National Laboratory. This study compares results in closely matched 
simulations from two staggered fermion actions, namely asqtad and p4fat3, and it is pro- 
viding the first large-scale simulation at N x = 8 with domain wall fermions. (See talks by 
M. Cheng [§, R. Gupta [§], and W. Soeldner [§].) 



• Chiral susceptibility. New insights into the behavior of the chiral susceptibility will change 
the determination of T c using this quantity [0]. 

• Equation of state. A new method has been proposed [j|]. 

• Transport coefficients. There are new ideas and methods for computing them [g]. 



• Spatial string tension. A new result agrees surprisingly well with 3D perturbation theory 

Ph. 



I will not have time to cover interesting studies of QCD-like theories with a large number of flavors 
[11]. And I regret that time and space did not permit covering all recent work in this field. 

I will try to give a general and fairly critical overview, using selected results from the parallel 
sessions as illustrations and leaving the details to the parallel sessions. After a brief review of 
lattice methodology I discuss potential cutoff problems with various actions and focus on issues 
and confusion in determining T c . Turning to results, I highlight some new methods and results for 
the equation of state, allude to recent progress in determining transport coefficients, and end with 
mention of a little surprise concerning predictions of dimensional reduction for the spatial string 
tension. 
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2. Lattice Methodology 

Lattice methods are especially well suited for simulating a quantum statistical ensemble in 
thermal equilibrium at fixed temperature T. We set a finite imaginary time interval 

aN r = l/T (2.1) 

for lattice spacing a and N t sites in imaginary time, and we impose periodic (antiperiodic) boundary 
conditions on the bosonic (fermionic) fields. Under these conditions the lattice Feynman path 
integral generates the quantum partition function for the underlying hamiltonian H, 

Z = Trexp(-///r), (2.2) 

in the continuum limit. Operator expectation values are thermal expectation values in this en- 
semble. Since simulations with standard methods are limited to equilibrium and near-equilibrium 
processes, to apply lattice results to the nonequilibrium conditions of heavy ion collisions requires 
phenomenological modeling. 

The temperature is varied by changing either N t or a. The latter strategy is more common. At 
fixed N t , decreasing the gauge coupling g 2 decreases a, so T grows. It is common now to adjust the 
bare lattice quark masses together with the lattice spacing so that zero temperature meson masses 
remain fixed. In this way variations in observables can be attributed to changes in temperature and 
not also to changes in the Hamiltonian. Such trajectories through parameter space are called "lines 
of constant physics." 

Of course, to connect with reality we need also to take the continuum limit. With the fixed 
N T strategy, the lattice is coarser at low temperatures and finer at high temperatures. For a given 
temperature, obviously, we approach the continuum by repeating the calculation at smaller a and 
larger N T . Contemporary lattice simulations have N T as large as 12 for some quantities [|i"2|], but 6 
and 8 are typical for expensive quantities such as the equation of state. Table [I] shows the relation- 
ship between lattice spacing and N r at T = 180 MeV, near the crossover temperature for QCD. We 
see that by standards of contemporary zero temperature simulations, thermodynamic simulations 
at N x = 6 and 8 are rather coarse at this temperature, and N T = 4 is extremely coarse. Thus we must 
be alert to the possibility of distortions due to cutoff effects. 





4 


6 


8 


10 


12 


a (fm) 


0.27 


0.18 


0.14 


0.11 


0.09 



Table 1: Lattice spacing vs. N r at T = 180 MeV. 

The continuum limit can be expensive. For the equation of state the computational cost grows 
with decreasing lattice spacing as a~ n . This places a high premium on reducing undesirable cutoff 
effects at a coarse lattice spacing. The degree of "improvement" of the lattice action is a significant 
factor. 

The most extensive recent simulations use staggered fermions with varying degrees of im- 
provement. The asqtad staggered formalism is designed to eliminate cutoff effects at &{a 2 ), leav- 



ing errors at G{a s a ) (see references in [13]). The p4fat3 staggered action is also improved, but 
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action 


A 2 /A 


A4/A0 


A 6 /A 


standard staggered 


248/147 


635/147 


3796/189 


Naik 





-1143/980 


-365/77 


p4 





-1143/980 


73/2079 


standard Wilson 


248/147 


635/147 


13351/8316 


hypercube 


-0.242381 


0.114366 


-0.0436614 


overlap/ 
domain wall 


248/147 


635/147 


3796/189 



Table 2: Continuum limit scaling behavior of free massless quarks in various lattice formulations, based 
on an expansion (2.3) of the pressure in powers of 1/N% from [[lq]. Shown are ratios of the expansion 
coefficients to the ideal, leading Stefan-Boltzmann coefficient. A small ratio indicates good scaling. 



it does not eliminate all such effects []14|]. Both actions also improve the free-quark dispersion re- 
lation, which is desirable in a high temperature deconfined environment. The Budapest- Wuppertal 
action with stout gauge links and unimproved staggered fermions does not improve the dispersion 



relation [15], but it does reduce taste-splitting effects. 

Improvement is good, but one may carry it too far. Improvement tends to fatten action oper- 
ators, in which case localization could become an issue. It is plausible that the extent to which a 



lattice simulation approximates the quantum partition function (2.2) depends on the locality of the 
lattice transfer matrix. Ideally the localization length of the action I should be much less than l/T 
or in lattice units, much less than N t . 

For free fermions cutoff effects for various lattice formulations can be studied analytically. 



Recently Hegde et al. [ |16| ] looked at deviations from the expected free-fermion Stefan-Boltzmann 
relation for the pressure p as a function of 1/A^ (equivalently a 2 ) and chemical potential fl/T: 

^=JLA 2k P 2k (^/KT)^-J (2.3) 

The leading term A is the Stefan-Boltzmann term. The ratios of higher coefficients A^IAq mea- 
sure the strength of the cutoff effects. These terms measure the ability of the action to approximate 
the continuum free fermion dispersion relation. Table ^ reproduces their results for a variety of 
actions. We see that the hypercube action JTv| ] has pleasingly small coefficients. The Naik (asqtad) 
and p4 (p4fat3) actions remove the second order term as designed, but the p4 action is better at 
sixth order. The standard (unimproved) staggered action (regardless of gauge-link smearing) does 
as poorly as the standard (and clover-improved) Wilson actions. The overlap and domain wall 
actions constructed from the standard Wilson kernel inherit its poor behavior. 

In selecting a fermion action for thermodynamics, there should be no excuse for deliberately 
building in poor continuum scaling. 

Another recent study confirms pronounced cutoff problems with free chiral fermion actions 
based on the standard Wilson kernel. Gavai and Sharma calculated the ratio of the lattice energy 



density to the expected Stefan-Boltzmann value for overlap and domain wall fermions [18]. Their 
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Figure 1: Deviation of the lattice free quark energy density from the Stefan-Boltzman continuum energy 
density as a function of N T from [^8J. (£ = N S /N T — 4). Left panel: overlap for various mass shifts M. Right 
panel: domain wall fermions for various L s at M = 1 .55. 




Figure 2: Splitting of the pion taste multiplet vs. a 2 a 2 for a wide range of lattice spacings. Splitting is mea- 
sured as the difference of the squared masses of the member and the Goldstone member. The plot symbols 
distinguish the members of the multiplet. The line is drawn with unit log-log slope to test proportionality to 
a 2 a 2 . 



results, reproduced in Fig. [I], show slow continuum scaling and an oscillation related to negative 
eigenvalues of the transfer matrix. 

Staggered fermions have the awkward problem of extra "taste" degrees of freedom. The stan- 
dard "fourth root" trick gives an approximately correct counting of species, but hadrons in the 
statistical ensemble still come in multiplets with a range of masses. This is especially visible in 
the pion spectrum. The splitting of taste multiplets is predicted to decrease as &(a 2 aj) in terms of 
lattice spacing and color fine structure constant a s . Figure || confirms this trend for the pion taste 
multiplet as a function of lattice spacing. 
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In judging how closely a simulation comes to the physical point, devotees of staggered fermions 
may be tempted to focus on the lowest member of the pion multiplet, i.e. , the Goldstone boson pion. 
While this practice is correct when the Goldstone boson can be an isolated external state, as in a 
zero temperature Green's function, in a thermal ensemble, all members of the multiplet participate. 
Thus it is more appropriate in thermodynamics simulations to compare an average multiplet mass, 
e.g. , the rms pion mass with the physical pion mass. The physical point is reached only by reducing 
the lattice spacing together with the light quark mass. It is incorrect to claim a thermodynamics 
calculation is done at a physical pion mass when the rms mass is still much higher. 

Is taste-symmetry breaking really a problem for thermodynamics? It is believed to be most 
dramatic for the pion and less noticeable for more massive states [19]. One could argue that close 
to the crossover temperature away from the critical point, so many excited states participate, as in 
the resonating hadron gas model, that pions do not matter much. But if approach the critical point 
at fixed lattice spacing, taste splitting is likely to have a strong effect on the critical behavior: we 
may even get a chiral symmetry restoring transition in the wrong universality class. And certainly 
at quite low temperatures where pions dominate the statistical ensemble, taste splitting makes a 
difference. 



3. Phase Diagram and Determination of T c 

3.1 Current consensus 

Is there a genuine phase transition separating a low temperature confined phase with spon- 
taneously broken chiral symmetry from a high temperature deconfined phase with restored chi- 
ral symmetry? The answer depends on the number of quark flavors and their masses. Figure || 
sketches the current qualitative theoretical consensus for the case of 2+1 flavors of quarks with 



masses m u = and m s [20]. In the upper right corner the quark masses are so large they play 
no role in the statistical ensemble, and we enter the well-studied regime of a first order confin- 
ing/deconfining transition in pure Yang-Mills theory. At low temperature chiral models predict a 
first order transition for degenerate masses, shown in the lower left cormer, and a second order tran- 
sition at large m s when m u = rrid = (provided the chiral anomaly does not vanish at the transition) 
dUl. 1 The low -mass first order region is bounded by a critical line above which the transition is 
only a crossover. 

Figure || is only qualitative. To say whether there is a phase transition at physical quark masses 
requires numerical simulation. The long-standing consensus has been that it is only a crossover. 



Aoki et al. have made a strong case for this conclusion [24]. Locating the actual critical line is 



challenging, since it occurs at small quark masses and is quite sensitive to cutoff effects [25, Eg]. 



DeForcrand and Philipsen ||27|, |28|] have done an impressively high statistics study that maps out 



the phase boundary, but only with unimproved staggered fermions at N T = 4. They have begun 



workatN T = 6 [28]. 



'in contradiction to this expectation, D'Elia, Di Giacomo, and Pica found indications of a first-order transition using 
an unimproved staggered fermion action and N T = 4 Jj^. It is important to check this conclusion with a more refined 
action. In support of this expectation, Kogut and Sinclair have found a second order phase transition, but in the 0(2) 
universality class, rather than 0(4) [^]. 
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Figure 3: Sketch from [|20p of the phase diagram for QCD at zero baryon density in 2 + 1 flavor QCD as a 
function of the light quark masses showing regions where a high temperature phase transition or crossover 
is expected. For a second-order phase transition, the universality class is shown. Whether the expected 
tricritical strange quark mass m 1 ™ is higher or lower than the physical strange quark mass mf 1 ^ is not yet 
firmly established. 

3.2 How precisely can we know T c l 

So what is the crossover temperature T c l One must ask, first, for what purpose do we need to 
know it? If there is only a crossover, the determination of T c is unquestionably imprecise, as the 



Budapest- Wuppertal group has emphasized [29]. For phenomenology it should be good enough 



to determine the temperature range over which an interesting quantity, such as the energy density 
changes rapidly. Each observable may give a somewhat different answer. 

For (presumably) unphysical quark masses for which a genuine phase transition occurs, the 
transition temperature T c is unambiguous and precision is achievable. The observables used to 
locate it must, obviously, have a sensible continuum limit, and they should expose the critical 
behavior. 

Finally, in determining the temperature there is a related question of setting the lattice scale. 
On coarse lattices, the result depends strongly on the physical quantity used to set the scale. 

With these preliminaries in mind we examine a variety of observables that have been used to 
determine the crossover temperature. First we consider "deconfmement-type" observables. 

3.3 Deconfinement observables 

3.3.1 Strange quark number susceptibility and equation of state 

The strange quark number susceptibility measures fluctuations in the strange quark number: 

% s = (Nf)/(VT) (3.1) 

The energy density and pressure are also good indicators of the progress of deconfinement. 
Recent results from the HotQCD collaboration are shown in Fig. ^ [|5, 3C]. 
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Figure 4: Left panel: Strange quark number susceptibility divided by the square of the temperature vs. tem- 
perature in MeV units (bottom scale) and ro units (top scale) foiN T = 6 and 8. Right panel: equation of state 
showing energy density and three times the pressure, both divided by the fourth power of the temperature vs. 
temperature for N T = 8. Measurements are taken along a line of constant physics with m u j = 0.1m s . These 
preliminary results are from a HotQCD study comparing p4fat3 and asqtad staggered fermion formulations 
[||30|. The blue error bars on the pressure curve indicate the size of the error. The magenta bar shows 
a systematic error from setting the lower limit of the pressure integration. The vertical bands here and in 
HotQCD figures below indicate a temperature range 185 - 195 MeV and serve to facilitate comparison. 



We note that the asqtad and p4fat3 results are in fair agreement. We see that there are dramatic 
changes in both observables over the temperature range 180 - 200 MeV. The scale has been set 
through the Sommer parameter ro (or r\), which, in turn, is calibrated by the measured splitting 
in bottomonium [31]. The rough agreement between the N t = 6 and 8 gives some support for the 
utility of this scale determination. 



3.3.2 Universal critical behavior 

To determine the crossover temperature more precisely, one may look for the inflection point 
in the quark number susceptibility or a peak in specific heat. As we have remarked, the result is 
unambiguous only for a genuine phase transition. Both of these quantities are derivatives of the 
free energy / = — TlogZ, which leads to a unified treatment of their critical behavior fl32"[]. At 



a critical point, the free energy can be decomposed into singular and analytic contributions. The 



singular part scales according to [ |32| , 



f s {T^ q )=b- l f s (tb 



,l/(2-«) 



t 2-a 



(3.2) 



where from charge symmetry and analyticity, the scaling variable depends on temperature and 
chemical potential through 



T-T c 



+ c 



(3.3) 



For 0(4) the critical exponent is a as —0.25. This relation can be used to predict singularities in 
quantities expressible as derivatives of the free energy. At jj. q = the light quark number suscepti- 
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bility is 

Wr 2 ~0^- a , (3.4) 



The singularity is weak and masked by analytic contributions. Its temperature derivative has a 
stronger singularity 

d ^ - r« (3.5) 



dT 

The same leading singularity is found in the specific heat 



d 2 f s 

C v ~ ~ t a (3.6) 



and in the quartic quark number susceptibility c\ = ((N*) — 3(N 2 )): 



d^fs 



4~^~t- a (3.7) 



3.3.3 Screening free energy of a static quark 



The Polyakov loop measures the screening free energy of a static quark [34]. Thus it is a 
phenomenologically interesting deconfinement-type observable. It is not the derivative of the free 
energy, so it might not reveal critical behavior. It is customary to renormalize it by removing a 



temperature-independent self-energy [35], leaving 



F q (T)=-T\og[P KDOrm (T)}. (3.8) 

This quantity is shown in Fig. |[ The related susceptibility (essentially the variance in this quantity) 
is sometimes used to locate the transition temperature, but in numerical simulation the peak in this 
susceptibility weakens with increasing N t . 

3.4 Chiral observables 

Next, we consider observables usually associated with a chiral phase transition. 

3.4.1 Chiral condensate 

The chiral condensate has an ultraviolet singularity for nonzero quark mass and (at least for 
Nf = 2) an infrared chiral singularity at zero quark mass m uc f. 



(\j/Y)(a,m ud ,T) ~ < 



ci/2(a,T)^/m^ + cim ud /a z + analytic T < T c 

c\m u d/a 2 + cgtnjf + analytic T = T c (3.9) 

c i m uc j /a 2 + analytic T > T c 



The ultraviolet singularity appears in perturbation theory at the one-quark-loop level. It is temper- 
ature independent. The infrared singularity occurs in the chirally broken phase T < T c [0]. It is 
seen in chiral perturbation theory at the one-pion-loop level. A square root is the thermal analog of 
a chiral log. At T c we have the expected critical behavior. 

The ultraviolet singularity in this quantity poses a problem for the continuum limit. Since 
the condensates for all flavors have this singularity, taking an the appropriate linear combination of 
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Figure 5: Renormalized screening free energy of a static quark (from the renormalized Polyakov loop) vs. 
temperature in MeV units (bottom scale) and ro units (top scale) for N T — 6 and 8 from a HotQCD study 
comparing p4fat3 and asqtad staggered fermion formulations [^J [30|] . Measurements are taken along a line 
of constant physics with m uc j = 0.1m s . 
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Figure 6: Chiral condensate difference ratio vs. temperature in MeV units (bottom scale) and ro units (top 
scale) for N T = 6 and 8 from a HotQCD study comparing p4fat3 and asqtad staggered fermion formulations 
[|J p0| ]. Measurements are taken along a line of constant physics with m uc j = 0.lm s . 



light and strange quark condensates removes it, and dividing by the zero-temperature value removes 
the multiplicative renormalization factor [32]: 



D l>s (T) = (ifryf)\ e -!^{ VV )\ s A o (T)=D ts (T)/D o (T = 0) 



(3.10) 



The difference ratio A^ s is shown in Fig. |6[ We see a dramatic drop in this quantity over approxi- 
mately the same temperature range over which we saw a rapid rise in Fig. ^ 
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3.4.2 Chiral susceptibility 

The chiral susceptibility measures fluctuations in the chiral condensate. Loosely speaking it is 

X = d{yf\i/)(a,m,T)/dm (3.11) 

where it is customary to distinguish the connected and disconnected contributions in terms of the 
Dirac matrix M: 



Xdisc — y 



T 

T 



' (TrAT 1 ) 2 ) - (TrAT u2 



Zconn = --(TrM- 2 ). (3.12) 

Tradition holds that a peak in #disc marks the crossover, but see the discussion of possible distortions 
below. 

Figure [7] shows recent results for the disconnected susceptibility. The peak occurs in roughly 
the same temperature range as the dramatic changes seen in the previous observables in Figs. ^| and 
|6| The peak height increases as expected as the light quark mass is decreased. 

The first exploratory N T = 8 results for the domain wall action [Q] are shown in Fig. [j| Because 
domain wall calculations are vastly more expensive, the domain wall effort is certainly not as 
advanced as other efforts. The L s = 96 study was undertaken to assure a small residual quark mass 
over the range of couplings shown. 



The derivative in Eq (3.1 1) generates the correlator of the relevant chiral condensates, inte- 



grated over the space-time volume: 



x = C(p = 0,T) = J d 4 xC{x,T) 
C(x,T) = (YW{x)WW(0)) (3.13) 



To be more precise, the derivative in the definition of the susceptibility (3.11) can involve any of 
the flavor condensates and any of the quark masses. For the light quark condensates, it is useful 
to distinguish the isosinglet and isotriplet chiral condensates, according to the isospin content of 
the operators in the correlator. These quantities are linear combinations of the disconnected and 
connected susceptibilities, namely, ^ sing = Xdkc + 2^ conn and Xtnp = 2Zconn- 

The ultraviolet and infrared singularities of the isosinglet chiral susceptibility can be easily 
inferred from Eq Q [0]: 



Zsing ~ < 



c\/a 2 + ci/2(a,r)/(2 v / ra) + analytic T < T c 

ci /a 2 + c s m 1 / s - l + analytic T = T C (3.14) 

c\/a 2 + analytic T > T c 



We see that it also suffers from an ultraviolet divergence, which makes it increasingly noisy as 
the lattice spacing is decreased. The infrared singularity at zero quark mass in the chirally broken 
phase (T < T c ) arises from the vanishing of the pion mass in that limit. For this reason at small 
quark mass we should expect not only a peak at the transition point, but we should expect the 
susceptibility to grow at lower temperatures as well 
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Figure 7: Left panel: Disconnected light quark susceptibility vs. temperature in MeV units (bottom scale). 
Right panel: closeup of the peak region. Lines merely connect the points. Red symbols, asqtad fermions. 
Blue symbols, p4fat3. Filled squares and circles are along a line of constant physics with m uc i = 0.1m s . Open 
circles, with m ud = Q.2m s , filled triangles with m ud — 0.05ms. All results are HotQCD preliminary [||6} |(]|. 




P 



Figure 8: Chiral condensate (left scale) and disconnected chiral susceptibility (right scale) vs. the gauge 
coupling parameter j3. Results are from a HotQCD study of 2 + 1 flavor domain wall thermodynamics at 
N T = 8 for two choices of L s ffl. Measurements are taken with light and strange quark masses fixed in lattice 
units (am m i = 0.003 and am s = 0.37). 
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Figure 9: Sketch of the expected behavior of the light screening spectrum for the %, fo and 2n threshold vs. 
temperature in units of the crossover temperature T c . 



To remove the ultraviolet singularity the Budapest/Wuppertal group subtracts the zero temper- 
ature value and multiplies by the square of the quark mass to cancel multiplicative renormalization 
factors: 

m 2 q \x(m g ,T)-x(m q ,0)]/T 4 (3.15) 

Compared with the uncorrected susceptibility, along a line of constant physics away from the chiral 
limit, this definition tends to shift a peak to lower T because m 2 /T 2 decreases with increasing T. 
If instead of a peak, there is a shoulder, it might induce a peak. 

3.4.3 Screening masses as indicators of the transition 



From Eq ( 3 . 1 3h we see that a spectral component of mass M{T) and weight p(M,T) in the 
correlator contributes p (Af , T) /M(T) 2 , which is singular when M(T) vanishes. Since the correlator 
is integrated over imaginary time, we can analyze the spatial dependence of the zero-Matsubara- 
frequency correlator to determine its spectrum. The masses in that case are "screening masses" 
[[36|]. Figure ^ sketches a possible scenario for the temperature dependence of the low spectral 
components of the isosinglet chiral condensate. In the chiral limit the fo must be degenerate with 
the pion for T > T c , so for a continuous transition, it must drop to zero there. As the light quark 
mass is decreased, the two-pion threshold also vanishes, leading to an infinite chiral susceptibility 



for T < T c as well, as reflected in Eq ( p. 14h . Close to the chiral limit, instead of a peak marking the 
transition, one might expect a cliff. 

At this conference on behalf of the RBC/Bielefeld collaboration Laermann reported new mea- 



surements of screening masses in the scalar and vector isotriplet channels []37[]. Some of their 
results are shown in Fig. [T(| We see that the % shows the behavior sketched in Fig. ^. Results for 
the /o are not available, but the ciq shows a steep dip at the crossover. 

As we see, screening masses may prove to be a useful indicator of the chiral transition, since 
they do not suffer from ultraviolet or chiral divergences, and they require no renormalization. 
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Figure 10: Screening masses for the pseudoscalar channel (upper panel) and scalar channel (lower panel) 
vs. temperature in a dynamical 2+1 flavor simulation with p4fat3 staggered fermions [ |38[ ] . Measurements 
were taken along lines of constant physics with m n ~ 220 MeV, mg = 500 MeV and N T = 6 and 8 [p7[]. 



3.5 Scale determination 

To quote the transition temperature T c in physical units requires a scale determination. The Bu- 
dapest/Wuppertal group favors setting the lattice scale with fy, whereas MILC and RBC/Bielefeld 
use the Sommer parameter rn or the related r\. At current typical lattice spacings and quark masses 
in staggered fermion simulations the fx scale results in a 10 to 20% lower temperature in MeV than 
the rn scale. This discrepancy vanishes at the physical quark mass and in the continuum. We should 
choose the scale so that the crossover temperature scales well. Of course, given the ambiguities 
in determining the crossover temperature, even in lattice units, that is an imprecise condition. The 
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deconfinement-type variables are more useful to phenomenology. For them the ro scale seems to 
give reasonably consistent results as the lattice spacing is decreased, as we can see from Figs. |] 



and g and Aoki et al. [ J29Q (Fig. 4). Thus there appears to be no reason to abandon the ro scale for 
now. 

In the past two years there have been some seemingly contradictory estimates of the transition 
temperature. Aoki et al. reported that at physical quark masses in the continuum, T c = 151 (3) (3) 
MeV from a peak in the chiral susceptibility and T c = 175(2) (4) from the inflection point in the 
quark number susceptibility and the Polyakov loop variable. These are to be compared with an 



older result from the MILC collaboration 169(12)(4) MeV based on the chiral susceptibility [39] 
and a more recent determination of 192(7) (4) MeV by the Bielefeld/RBC group based on a com- 
bination of chiral and deconfining observables pOp. The last two groups used the ro (or r\) scale. 
The Budapest/Wuppertal group has carefully listed sources of the discrepancy, which include am- 
biguities in locating the crossover, their preferred renormalization of the chiral susceptibility, and 
their preference for the fg scale. This year we can add to the list the possibility that the chiral 
susceptibility develops an asymmetric peak or shoulder, which would be even more sensitive to 
the renormalization procedure and should not be modeled by a parabola. A combination of these 
effects could certainly account for the discrepancy. 

4. Equation of state 

The equation of state is fundamental to hydrodynamic calculations of the expansion of hot 
hadronic matter. 

4.1 Standard integral method 

The currently popular method for calculating the equation of state begins with the lattice- 
thermodynamic identity at fixed N t , which expresses the trace of the energy momentum tensor or 
"interaction measure" / in terms of the derivative of the log of the partition function Z: 

TdlnZ 

I = e -3p = --——. (4.1) 
V ama 

The derivative with respect to lattice spacing is taken with fixed output hadron masses. Thus it 
involves the derivative of the bare lattice parameters, i.e., the gauge coupling g, quark masses, and 
for some actions the tadpole coefficient uq, with respect to the cutoff scale and the expectation 
values of the action operators. All of these nonperturbative quantities are readily calculated in 
lattice simulations. 

We normalize the energy and pressure to zero at zero temperature. This eliminates an ultravi- 
olet divergence of order I /a 4 . This is done by subtracting the zero temperature quantity, calculated 
with the same bare parameters. Because we are subtracting two ultraviolet singular quantities, we 
must increase the simulation sample size dramatically as the lattice spacing is decreased. The cal- 
culational cost thus grows steeply as we approach the continuum limit. In what follows we will 
assume this subtraction has been done for all thermodynamic quantities. 

A second identity becomes 

p dlnZ 



T dV 



(lnZ)/V (4.2) 
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Figure 11: Details of the dependence of the interaction measure on temperature in MeV units (bottom scale) 
and ro units (top scale) for three temperature ranges left to right: low, mid, and high, for N T — 6 and 8 from 
a HotQCD study comparing p4fat3 and asqtad staggered fermion formulations [M 0, 63]. Measurements in 
most cases are taken along a line of constant physics with m u( ( = Q.\m s . Results in the high temperature 
range at N T — 4 are from [5^]. In the low temperature range the magenta curve is the prediction of a hadron 
resonance gas model. The other curves in that range are spline fits to the data. The curves in the high 
temperature range are fits to a quadratic in 1/T 2 . 



in the thermodynamic limit for which InZ oc V. Finite-size deviations from this limit could produce 



deviations from the Stefan-Boltzmann law p41|, |42|]. Such finite size effects may be important at 
ultrahigh T, where we would like to compare with perturbation theory p3|]. 



Putting Eqs (4.1) and (4.2) together gives the the pressure as the integral of the interaction 
measure from coarse to fine lattice spacing, i.e., low to high temperature: 



Vp 
T 



Vp 




Ana 


~ ~T 

a 




JXnan 



Ana y> 

/ -(e'-3pVW. 

Jinan ■* 



(4.3) 



If the lower limit is sufficiently low in temperature, the pressure is zero. With pressure and interac- 
tion measure in hand we immediately get the energy density and entropy. 

Recent results for the interaction measure for N t = 6 and 8 are shown in Fig. 11. The resulting 
equation of state and pressure were shown in Fig. |] and the entropy density is shown in Fig. 12. 
Whether there is a statistically significant disagreement between the asqtad and p4fat3 results in the 
central region remains to be determined after further data are accumulated. In the high temperature 
range T G [250,700] MeV, the results can be fit to 



(e-3p)/T 4 = b/T 2 + c/T 4 



(4.4) 



There seems to be no need yet to include perturbative 1 / log T terms in the fit coefficients b and 
c. The plot in the low temperature range compares the lattice result with predictions of a hadron 
resonance gas model. Since the lattice calculations include cutoff effects, it would be premature to 
draw conclusions based on a disagreement at this level. 

4.2 New T integral method 

In the past year the WHOT collaboration introduced a new method in which the integral over 
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Figure 12: Entropy density divided by the third power of the temperature vs. temperature in MeV units 
(bottom scale) and ro units (top scale), for N T = 6 and 8 from a HotQCD study comparing p4fat3 and asqtad 
staggered fermion formulations [|^, [30| ], Measurements in most cases are taken along a line of constant 
physics with m uc / = O.lrn,. 



lattice spacing at fixed N t is replaced by an approximate integral over 7V T at fixed lattice spacing 
(i.e., fixed bare lattice parameters) [||, |4|]. 

The method starts from an alternative form of Eq ( |4.1[ ) with Eq (1.2): 



I/T 4 = d ~^. (4-5) 

The pressure is then computed by integrating the interaction measure with respect to In T or 
equivalently ln(N t ). Since N z is an integer, to reduce discretization errors in the sampling of the 
integrand, one must reduce the temporal lattice spacing a t . An anisotropic lattice helps. 

With this potentially computationally cheaper method the zero temperature subtraction is com- 
mon to all N x , and with bare lattice parameters fixed, one necessarily works along lines of constant 



physics. Figure 13 shows the result of a test calculation done for pure SU(3) Yang-Mills theory. 



5. Plasma Structure 



In addition to the phase diagram and equation of state, lattice simulations provide information 
about the structure of hot hadronic matter. I mention two recent developments: a new effort to 
determine the shear and bulk viscosity over a range of temperatures and a new study of the spatial 
string tension. 

5.1 Transport Coefficients 

Analysis of RHIC heavy ion collisions suggests that high temperature hadronic matter is an 
exceptionally good fluid. To confirm this hydrodynamical characterization requires computing 
the transport coefficients, namely the shear (tj) and bulk (£) viscosities. They are obtained from 
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Figure 13: Equation of state (interaction measure, energy density and pressure) for pure Yang-Mills theory, 
obtained using the T integral method at fixed lattice spacing a a = 0.097 fm and aspect ratio a T /a a = 4 
*4l 



correlators of the energy-momentum tensor at temperature T 

C(x ,x,T) = (T^ v (x Q ,x)T pa (0).) (5.1) 
We need its spectral function p, which we obtain from the Kubo formula 

I"" , s cosh co(x - 1/27) 

The transport coefficients are obtained from the low-frequency behavior of the spectral function 
nm -*^ e!£2!gm C(r) = | lim ^M) (5 .3, 

m^o CO 9 ffl^O CO 

This has been a well known challenging problem since it was first attempted by Karsch and Wyld 
[E5l]. The correlator is noisy, requiring high statistics. Going from a Euclidean correlator C(xo) 
to p(co) is a very difficult inverse problem. Because of time-reflection symmetry, a simulation at 
N t = S has only five, typically noisy, independent values. 

Possible remedies include (1) assuming a functional form for p and fitting its parameters, (2) 
decreasing the time interval a t , allowing a larger 7V T , and (3) adding further constraints on p, such 
as maximum entropy. 

Meyer [^, 4^, 47] has done a new high statistics calculation in pure Yang -Mills theory and uses 



a paramerization of the spectral function in terms of an optimized basis set that folds in appropriate 
perturbative behavior. For the ratio of shear viscosity to entropy density, he finds f]/s = 0.134(33) 
at 1.65r c where perturbation theory gives 0.8, and for the ratio of bulk viscosity to entropy density, 
(,js < 0.15 at 1.6571- and £/j < 0.015 at 3.2r c . 

5.2 Spatial string tension 

Despite its popular characterization as deconfined, high temperature hadronic matter retains 
vestiges of confinement. Space-like Wilson loops still exhibit the area-law behavior associated 
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Figure 14: Temperature divided by the square root of the spatial string tension o s vs. temperature in units 
of the crossover temperature 7b (lower scale) and in ro units (upper scale) for 2+1 flavors of p4fat3 quarks 
on lattices with N T — 4, 6 and 8. The solid curve (with uncertainties indicated by the dashed lines) is the 
prediction of the dimensionally reduced theory Jl0| ] 



with confinement. This is readily seen by considering dimensional reduction, in which for T S> T c 
the short Euclidean time dimension is collapsed, leaving three spatial dimensions, one of which is 
reinterpreted as the Euclidean time coordinate of a 2+1 dimensional field theory. 
The reduction of QCD has these characteristics: 



Quarks acquire a large 3D mass . / (tiT) 2 + m- 



9 

• The fourth component of the color vector potential Aq becomes a scalar field, and we get a 
confining gauge-Higgs theory. 



The 3D and 4D gauge couplings are related through g3 = gA^/T. 
i The spatial Wilson loop gives the 3D potential and 3D string tension. 
In a recent calculation Cheng et al compared the spatial string tension of the full 4D theory 



with its predicted behavior in 3D perturbation theory QIQ]. The comparison is shown in Fig. [14. 
The good agreement with perturbation theory at temperatures as low as \ .5T C is unexpected. 



6. Conclusions 

In a reasonably well-matched calculation, new high statistics results from HotQCD show good 
agreement between two different staggered fermion formulations, i.e., p4fat3 and asqtad. Not 
surprisingly, simulations with these inexpensive algorithms are more advanced than those with 
other fermion actions, as we have seen from the first exploratory domain-wall-fermion simulations 
at N t = 8 with a quite small residual quark mass. To make progress we need to understand the 
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importance of cutoff effects and to come closer to the physical point. Calculations with other 
fermion actions can provide an important check, but those actions must be improved at least to the 
same level as the staggered fermion actions before they can play this role effectively. 

We are learning more about the phase structure of zero-baryon-density QCD as a function of 
the light quark masses, but these results are especially sensitive to cutoff effects. More work is still 
needed. 

There has been recent progress in methodology. The WHOT collaboration has developed 
a new method for determining the equation of state, and Meyer has proposed new methods for 
determining transport coefficients. 

Finally, in measurements of the spatial string tension, we have seen interesting agreement with 
predictions of dimensional reduction. 
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